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Abstract 

Relative to western populations, the percentage of women diagnosed with breast cancer at a young age in Lebanon is high. 
While the younger age of the Lebanese population compared to the West certainly contributes to this difference, potential 
genetic, reproductive and/or biological factors likely play an important role. The objective of this study is to investigate the 
contribution of miRNAs in this setting through the analysis of the expression of five reported dysregulated miRNAs, miR- 
148b, miR-lOb, miR-21, miR-221, and miR-155 in 20 normal and 57 cancerous breast tissues from Lebanese breast cancer 
patients. After finding their relative expression by quantitative reverse transcription real time PCR, the results were analyzed 
with respect to the patients' clinical and histopathology presentations. Compared to normal breast tissues, significant 
upregulation of miR-155, miR-21 and miR-148b, notable downregulation of miR-lOb and non-significant expression of miR- 
221 were observed in tumor tissues. Moreover, miR-lOb was significantly underexpressed in estrogen/progesterone 
receptor (ER/PR) negative tumors relative to ER/PR positive tumor tissues. miR-155 was also significantly overexpressed in 
postmenopausal patients and in those of age at diagnosis greater than 40 years old as well as in PR negative or in human 
epidermal growth factor 2 (Her2) positive tissues. This study is the first one to report miRNA expression patterns in Lebanese 
breast cancer patients. We found that differential miRNA expression in breast cancer could be variable between Lebanese 
and Western populations. miR-lOb was positively correlated with the ER and PR status and miR-155 could be a noteworthy 
biomarker for the menopausal state, age at diagnosis, PR and Her2 status. Hence, miRNA can be used as biomarkers for early 
breast cancer detection. 
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Introduction 

Breast cancer is one of the leading health concerns worldwide, 
affecting over one million women every year. In Lebanon, it is one 
of the most common type of cancer constituting about one third of 
all female cancers [1]. Interestingly, a significant number of 
Lebanese breast cancer patients were noted to be of young age at 
the time of diagnosis as 22% of the cases were below the age of 40 
years old compared to 6% in the Western populations [2]. 
Moreover, Lebanese women who are diagnosed at young age (less 
than 35 years old) and in their premenopausal state were shown to 
present with a more aggressive disease and poorer survival in spite 
of adequate therapy [3] . The presence of signs of more aggressive 
features in breast cancer in young women, and the occurrence of 
breast cancer in young Lebanese women 10 years earlier in age 
than those in the West, strengthens the importance of determining 



the biological factors behind those differences, and perhaps 
revealing novel biomarkers for early screening and detection of 
breast cancer. 

Since the discovery of microRNA (miRNA) in C. elegans twenty 
years ago, this major subclass of non-coding RNA molecules act as 
gene modulators mostly at the posttranscriptional level by causing 
translation repression or degradation of mRNA [4] . miRNAs play 
diverse roles in normal cellular processes such as cell cycle, 
proliferation and apoptosis as well as in disease conditions 
including cancer, diabetes, neuro-degenerative disorder and 
cardiac hypertrophy [5,6] . miRNA was first correlated with breast 
cancer by Iorio and his colleagues. Using microarray analysis, they 
discovered a differential miRNA profile pattern between cancer- 
ous and normal breast tissues [7]. miRNAs were later shown to 
modulate tumor suppressor and oncogenic pathways thereby 
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contributing to the different stages of breast cancer and acting as 
regulators of cell cycle progression, apoptosis, angiogenesis, 
epithelial to mesenchymal transition, tumor microenvironment, 
migration, invasion and metastasis [8] . Previous studies attempted 
to correlate the different miRNA expression profiles in tumor 
tissues with the histological diagnosis stage and other variables [9] . 
Interestingly, studies have also shown that miRNAs are not only 
present in tissues, but also in biological fluids including plasma. 
This suggests that miRNAs have the potential to serve as non- 
invasive diagnostic and prognostic biomarkers as well as possible 
therapeutic targets. 

The oncogenic role of miR-155, miR-21, miR-lOb and miR- 
221 and their potential use as biomarkers in breast cancer have 
been reported extensively [10-13]. In addition, circulating miR- 
148b was found to have a discriminatory potential for young 
breast cancer patients [14]. This study aims to investigate the 
expression patterns of these miRNAs, which are known to be 
dysregulated in breast cancer, in Lebanese breast cancer patients 
of varying clinical and histopathology presentations. 

Materials and Methods 

Breast cancer tissue specimens 

After obtaining the approval of the Institutional Review Board 
(IRB) at the American University of Beirut, formalin fixed paraffin 
embedded (FFPE) sections from invasive ductal carcinoma 
specimens (N = 57) and normal adjacent tissues (NAT) (N = 20) 
were obtained at the American University of Beirut Medical 
Center (AUBMC) in Lebanon from samples that were collected 
between 1997 and 2012 and were previously recruited for a breast 
cancer study [15]. Patients who participated in the originally IRB 
approved study signed a written informed consent form. Those 
who agreed that their stored breast cancer FFPE tissues, if 
available, are used for further research were included in this study. 
The Pathology department at AUBMC provided four to five 
20 nm ribbons from the tumor biopsy core and five 20 um ribbons 
from the adjacent normal tissue when available. Sections were 
taken from paraffin embedded tissue blocks containing at least 
90% viable invasive tumor cells. The corresponding adjacent 
normal breast tissue sections included epithelial mammary 
parenchyma and were entirely free of significant pathology such 
as atypical hyperplasia or carcinoma. Benign proliferative changes 
such as adenosis were allowed. Clinical and pathological data, 
such as the age at diagnosis, menopausal status, tumor pathology, 
stage, grade, ER status, PR status, and Her2 overexpression was 
available for the 57 tumor samples. 

Total RNA extraction 

Total RNA was extracted from 20 normal and 57 tumor tissues 
in accordance to the protocol of RecoverAll Total Nucleic Acid 
Isolation Kit for FFPE samples (Ambion, USA). Briefly, FFPE 
samples were first deparaflinized using xylene at 50°C and then 
washed twice with ethanol to remove the xylene. Proteins were 
digested by incubation of the samples with protease enzyme for 
15 min at 50°C and then at 80°C. Total RNA was then captured 
in glass-fiber filter cartridges with the help of an isolation additive 
mixture followed by washing steps with high ethanol-wash buffers. 
After DNA digestion using DNase for 30 min, RNA was purified 
by washing and elution. RNA concentration and quality was 
assessed using the Nanodrop ND1000 and then stored at — 80°C. 
Only high quality samples were used for downstream applications. 



miRNA expression by quantitative real time-polymerase 
chain reaction (RT-qPCR) 

Reverse transcription of ten nanograms of the total RNA was 
performed using the TaqMan MicroRNA Reverse Transcription 
Kit (Applied Biosystems, USA) according to the manufacturer's 
instructions. Small nuclear RNA RNU6B, human hsa-miR-16, 
hsa-miR-21, hsa-miR-221, hsa-miR-148b, hsa-miR-lOb, and hsa- 
miR-155 primers and probes were purchased as part of the 
TaqMan microRNA Assays Kit (Applied Biosystems, USA) with 
validated efficiency. cDNA synthesis was carried out in a multiplex 
reaction set up whereby two miRNA primers (for example miR- 
148b and miR-lOb) were used in each reaction with the 
endogenous control, miR-16 or RNU6B. RT-qPCR was per- 
formed using BioRad CFX96 Real Time System, CI 000 Thermal 
Cycler (Germany). Reactions using 10 ul of 2x TaqMan Universal 
Master Mix with no Amperase Uracil N-glycosylase (UNG) 
(Applied Biosystems, USA), 1 ul of the corresponding 20x 
microRNA probe, 4 ul of DEPC treated water, and 5 ul of 
cDNA were performed in duplicates for each miRNA probe. 
cDNA Synthesis and RT-qPCR were repeated twice for each 
sample and each plate included: no reverse transcription control 
(NRT), no template control (NTC) and normal breast tissue 
samples. The normalization of the tumor tissues was based on the 
normal tissues present in the RT-qPCR plate to ensure inter-run 
calibration. The cycling conditions were 95°C for 10 min and 40 
cycles of 95°C for 15 seconds and an annealing temperature of 
60°C for 60 seconds. Using the AACt equation, the relative 
expression of the experimental miRNA was determined in the 
tumor samples compared to the normal tissue specimens using 
miR-16 or RNU6B as an endogenous control. 

Statistical Analysis 

Statistical Analysis was performed using SPSS software package 
version 18. Wilcoxon's rank sum test for one sample was used to 
compare the miRNA expression in the tumor versus the normal 
tissues. Mann-Whitney V nonparametric test was used for 
comparing two different groups such as premenopausal and 
postmenopausal, but if more than two groups were compared, 
Kruskal- Wallis nonparametric test was used. Spearman's rho (r) 
was calculated to find a correlation between two variables. A p- 
value<0.05 was considered statistically significant. 

Results 

Baseline Demographics 

The clinical and pathological data have been retrieved for all 57 
samples of the breast cancer patients (Table 1). The breast cancer 
pathology was invasive ductal carcinoma for all the cases with 
78.9% ER positive, 66.7% PR positive and 33.3% Her2 
overexpression (defined as score 3 by IHC or score 2 by IHC 
with positive FISH) [16]. Most of the tumors were moderately to 
highly differentiated (42.1% and 47.4% for tumor grade 2 and 3). 
52.6% of the tumors were T2 with local growth and an average 
tumor size between two to five centimeters, 53.8% had lymph 
node involvement, and most had no distant metastasis (M0: 93%) 
according to the American Joint Committee on Cancer's (AJCC) 
Staging [17]. 

miRNA expression in breast cancer and NATs 

In order to compare normal and tumor breast tissues in the 
expression of the selected miRNAs (miR- 1 48b, miR- 10b, miR-2 1 , 
miR-221, and miR-155), reverse transcription real time PCR was 
performed on cDNA transcribed from RNA that was extracted 
from 20 normal and 57 tumor breast formalin fixed paraffin 
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embedded tissues. After normalization in the RT-qPCR with miR- 
16 or RNU6B as endogenous control, Wilcoxon's rank sum test 
was used to identify the miRNAs that were significandy different in 
the normal and tumor samples. Similar results were obtained upon 
normalizing to either one of the endogenous controls. miR-148b, 
miR-21 and miR-155 were significandy overexpressed (p<0.05) in 
the tumor tissues versus the normal tissues. While miR-221 
expression was not significantly different, miR-lOb was found to 
be significantly underexpressed in tumor as compared to normal 
tissues (Figure 1). Table 2 also shows the expression of the five 
miRNA in tumor versus normal breast tissues in sample groups of 
specific menopausal status (premenopausal or postmenopausal) or 
of specific age at diagnosis (younger than 40 or older than 40 years 
old). The expression results of these miRNAs in these groups were 
the same as when taking the total sample size with miR-148b, 
miR-21 and miR-155 significantly overexpressed and miR-lOb 
significantly underexpressed in tumor as compared to normal 
tissues (Table 2). 



miRNA expression and menopausal state or age at 
diagnosis 

The variation in the miRNA expression according to meno- 
pausal status and age group (£40 years old and >40 years) was 
analyzed by Mann-Whitney U test. When looking at the 
menopausal status or age at diagnosis, no significant difference 
in the miRNA expression was observed between the premeno- 
pausal and postmenopausal groups or between the age groups for 
miR-21, miR-lOb, miR-148 and miR-221. miR-155 was signif- 
icantly upregulated in postmenopausal patients when compared to 
premenopausal patients, as well as in the age group greater than 
40 years as compared to that of less than or equal to 40 years old 
(Figure 2). 

miRNA expression and clinical histopathological features 

To evaluate whether certain miRNAs are dysregulated with the 
various clinical and histopathological features associated with 
breast cancer, we analyzed our data by nonparametric statistical 
tests (Mann-Whitney U and Kruskal-Wallis tests) after dividing the 
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Figure 1. miRNA expression is dysregulated in breast cancer tissues in the Lebanese population. This boxplot displays the fold change 
of relative expression of miR-148b, miR-10b, miR-21, miR-221 and miR-155 in 57 tumor versus 20 normal adjacent breast cancer tissues that is 
measured by RT-qPCR with miR-1 6 used as an endogenous control. Whiskers represent minimum and maximum, the top, the bottom and the band in 
the box represent the first and third quartile and the median respectively. *denotes p<0.05, using Wilcoxon's signed-rank test. 
doi:1 0.1 371 /journal.pone.01 07566.g001 
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samples into groups with differential ER or PR expression or Her2 
overexpression, lymph node involvement, tumor size and grade. 
No analysis was performed on the different histological types since 
all tumor samples were of ductal carcinoma type, the most 
common type in the Lebanese population. Upon ER grouping, 
there were significantly lower miR-10 expression levels in ER 
negative specimens relative to the positive ones (Figure 3a). Upon 
PR grouping, data revealed a significant upregulation of miR-155 
and downregulation of miR- 1 Ob in the PR negative samples when 
compared to the PR positive samples (Figure 3b and c). As for 
Her2 overexpression grouping, only miR-155 was significantly 
upregulated in the Her2 overexpressed samples as compared to 
those not overexpressing Her2 (Figure 3d). Spearman-Rho test 
was out carried out to find correlations between the miRNA 
relative expression and the different clinical histopathological 
features. miR- 10b was positively correlated with the ER and PR 
status with correlation coefficients of r — 0.368 and r = 0.262 (p< 
0.05 for each) respectively and miR-155 was positively correlated 
with Her2 expression and negatively correlated with PR expres- 
sion with correlation coefficients of r — 0.308 and r=— 0.3 (p< 
0.05 for each)respectively. Lymph node involvement, tumor size 
and grade groups did not reveal any significant difference in the 
expression of the selected experimental miRNAs (Table 1). 

Discussion 

Compared to western populations, Lebanese women are more 
likely to be diagnosed with breast cancer at a younger age, which 
has been reported to be associated with adverse pathological 
features such as histological high grade, negative hormone 
receptors, HER2 overexpression and advanced stage at presenta- 
tion and hence carries a poorer prognosis [18]. Recent research 
explaining breast cancer at young age is focusing on possible 
molecular aberrations such as RANKL, PI3K pathways in order 
to identify possible therapeutic targets [19,20]. miRNAs are now 
widely investigated in different types of cancers and particularly in 
breast cancer. Their differential expression may be useful for the 
study of etiology of breast cancer, prognosis, and possible 
therapeutic targeting. This study looks at dysregulated miRNA 
profiles for Lebanese breast cancer patients where a significant 
number of patients are of young age. In this work, we have studied 
the expression of specific miRNAs (miR- 148b, miR- 10b, miR-21, 
miR-221 and miR-155) by RT-qPCR in cancerous and normal 
adjacent breast tissues and investigated whether their expression is 
linked to this increasing frequency for early onset of breast cancer 
in the Lebanese women. Our data showed that miR- 148b, miR-21 
and miR-155 were significantly overexpressed in tumor tissues 
versus normal tissues. While miR-221 expression was not 
significantly different, miR- 10b was also found to have signifi- 
cantly decreased expression in tumor compared to normal tissues 
(Table 2). The upregulation of miR- 21 and miR- 155 in tumor 
versus normal tissues was consistent with the literature and their 
S — g oncogenic role [7] . miR-2 1 has been reported as an oncomir that 

jj < ? contributes to motility invasion and metastasis by targeting tumor 

c ^ 'I * suppressors such as phosphatase and tensin homolog (PTEN), 

•| S a; 2 tropomyosin 1 , programmed cell death 4 (PDCD4) and mammary 

s £ t serine protease inhibitor (Maspin) [21-23]. miR-155 has been also 

S | o q. shown to be an oncomir acting as apoptosis suppressor and 

p E inducer of cell proliferation and chemoresistance through targeting 

~v 2 2 .o Caspase3, suppressor of cytokine signaling 1 (SOCS1) and 
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| o .s 5 [24,25]. On the contrary, the upregulation of miR- 148b was not 

E § | | in accordance with its reported function as a tumor suppressor that 

was downregulated in breast tumors [26]. However, our results are 



5 



September 2014 | Volume 9 | Issue 9 | e1 07566 



miRNA in Young and Lebanese Breast Cancer Women 



a) 



x 10- 



miR-155 




b) 



5 15 



miR-155 



oo° 



Figure 2. miR-155 is upregulated in breast cancer tissues of postmenopausal women and those greater than 40 years old. This 
dotplot shows the fold change of the significantly differentially expressed miRNAs upon dividing patients into groups based on the age at diagnosis 
and menopausal status. Relative expression of miR-1 55 in breast cancer tissues from patients of a) age at diagnosis less than or equal to 40 years and 
age greater than 40 years or of b) premenopausal and postmenopausal status was measured by RT-qPCR with miR-1 6 used as an endogenous 
control. The plot represents the mean with the standard error of mean as error bars. *denotes p<0.05 for tumor versus normal using Wilcoxon's 
signed-rank test. #denotes p<0.05, using Mann-Whitney U test. 
doi:1 0.1 371 /journal.pone.01 07566.g002 
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Figure 3. Relation of miR-155 and miR-1 Ob with the receptor profile. This dotplot shows the fold change of the significantly differentially 
expressed miRNAs upon dividing patients into groups based on ER, PR or Her2 status. Relative expression of a) miR-1 0b in ER negative and positive 
breast cancer tissues b) and c) miR-1 0b and miR-155 in PR negative and positive breast cancer tissues respectively d) miR-155 in Her2 negative and 
positive breast cancer tissues was measured by RT-qPCR with miR-1 6 used as an endogenous control. The plot represents the mean with the 
standard error of mean as error bars. *denotes p<0.05 for tumor versus normal using Wilcoxon's signed-rank test. #denotes p<0.05, using Mann- 
Whitney U test. 
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in agreement with the increased expression of miR- 1 48b obtained 
in ovarian carcinoma [27]. Interestingly, miR- 148b was also 
shown to be upregulated in the serum of young-aged breast cancer 
patients [14]. Downregulation of miR- 10b, an oncomir with 
metastatic potential, was obtained in metastasis free patients [7,28] 
and that is the case with patients in this study who mostly have no 
metastasis. As for miR-221, no significant difference in its 
expression was obtained in Lebanese breast cancer patient 
population which contrasts previous reports that demonstrated 
its increase in breast cancer due to its oncogenic role by regulating 
proliferation, angiogenesis and multiple tumor suppressors 
[10,29]. Hence, differentially expressed miRNA between normal 
and tumor breast tissues could be variable between Lebanese and 
Western populations. Remarkably, it was previously found that 
there is little overlap between the differentially expressed miRNAs 
identified from Caucasian American breast cancer patients and 
those from African American groups, suggesting a potential ethnic 
difference in miRNAs [30]. 

The significant downregulation of miR- 10b in ER/PR negative 
samples as compared to ER/PR positive samples and its positive 
correlation with ER and PR status establishes a relationship 
between miR- 10b and ER and PR expression. In addition, the 
significant upregulation of miR- 155 in PR negative samples 
compared to PR positive and its negative correlation with PR 
status were in agreement with a published study by Lu et al. in 
2012. Similarly, the significant overexpression of miR-155 in Her2 
positive relative to Her2 negative samples and its positive 
correlation with Her2 status were similar to what was previously 
reported [31,32]. miR-155 could not only predict the PR and the 
Her2 status but also its significant upregulation in postmenopausal 
samples and in patients of age greater than 40 years old made it a 
considerable breast cancer biomarker for postmenopausal patients 
or for those greater than 40 years old. On the other hand, finding 
differentially expressed miRNA specific for premenopausal 
patients or for those less than or equal to 40 years would be 
more beneficial for diagnosis of early onset breast cancer in 
Lebanese patients. Even though many studies have correlated the 
miR-2 1 expression with lymph node involvement and histological 
grade, our study did not show any statistically significant increase 
that could be due to the small sample size [33,34]. 

It is important to note that the NAT controls used in this study 
were not accurate normal breast tissues from healthy subjects but 
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microarray analysis to identify differentially expressed miRNA in 
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In conclusion, this study is the first one to report miRNA 
expression profiles in a Lebanese breast cancer population, a 
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miR- 10b was positively correlated with the ER and PR status. 
Moreover, miR-155 was shown as a noteworthy biomarker for the 
menopausal state, PR and Her2 status (Table 1). This suggests that 
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